Bose Hubbard model in the presence of Ohmic dissipation 
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We study the zero temperature mean-field phase diagram of the Bose-Hubbard model in the 
presence of local coupling between the bosons and an external bath. We consider a coupling that 
conserves the on-site occupation number, preserving the robustness of the Mott and superfluid 
phases. We show that the coupling to the bath renormalizes the chemical potential and the in- 
teraction between the bosons and reduces the size of the superfluid regions between the insulating 
lobes. For strong enough coupling, a finite value of hopping is required to obtain superfluidity 
around the degeneracy points where Mott phases with different occupation numbers coexist. We 
discuss the role that such a bath coupling may play in experiments that probe the formation of the 
insulator-superfluid shell structure in systems of trapped atoms. 

PACS numbers: 05.30 Jp, 03.75Lm, 03.75Hh 
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The Bose-Hubbard model (BHM) describes a wide va- 
riety of physical situations for appropriate parameter 
ranges Josephson junction arrays, ultracold atoms 
in optical lattices @, y], and polaritons or photons in 
arrays of resonant optical cavities 0, A key fea- 
ture of the BHM is the existence of a transition between 
Mott insulator and superfluid phases, recently observed 
for cold atoms in optical lattices This superfluid- 

insulator transition (SIT) is well understood theoretically 
@; E3, 11 1 and occurs as the magnitude of the inter-site 
hopping J is varied relative to the on-site repulsion U . At 
T = 0, there are insulating Mott phases for small enough 
J/U which exist in well-defined domains (Mott lobes) in 
the zJ/U-[j,/U plane (/x is the chemical potential and z 
is the co-ordination) , while in the large J limit there is a 
superfluid state with fully delocalizcd particles. 

In some, or all, of the various systems that are de- 
scribed by the BHM, it is likely that there will be cou- 
pling to external degrees of freedom in the environment 
that can be modelled as a bath that lead to equilibration 
in the BH system. At small but finite temperatures, the 
insulating nature of the Mott phase does not change, al- 
though the on-site occupation numbers undergo thermal 
fluctuations leading to finite compressibility [13]. Ex- 
periments on optical lattice systems are in this regime, 
where temperature is important (l3| . However, it is still 
a matter of debate whether heating takes plac e as the 
depth of the lattice potential is increased FJ, |l5|. It 
has been suggested that the main cause for such a tem- 
perature rise is the increase of the excitation energy gap 
combined with the necessity of maintaining a constant 
entropy [HI, 3 12 1- If heating can be important, then 
maintaining the system at a constant temperature may 
require the presence of some source of dissipation to re- 
move energy. The issue of dissipative dynamics in cold 
atom systems has attracted considerable recent atten- 
tion, in the contexts of cooling atoms in the lowest Bloch 
band 18], Bose- Fermi mixtures [l9| . the decay of su- 



percurrents 

Sip, and one dimensional cold molecular 
gases [13]. Hence, it is of great interest to address theo- 
retically the issue of how coupling to a heat bath affects 
the phase diagram of the BHM. 

In this Letter we investigate the influence of coupling 
to a heat bath (assumed to be Ohmic) on the standard 
phase diagram of the BHM [3, [Hi EH . We assume that 
the system of bosons is in full equilibrium with the bath 
at fixed temperature. Throughout our treatment, we 
assume for simplicity that T = 0, but briefly mention 
possible changes in the phase diagram due to finite tem- 
perature. We consider a coupling between the bath and 
bosons through the particle operator hi, implying that 
the number of bosons in the system is fixed. Hence the 
Mott insulator may still be defined as a state with in- 
teger expectation value for the occupation number no 
on each site, and in the superfluid phase there is a fi- 
nite phase stiffness. The coupling to a bath leads to: i) 
Renormalization of /i and U to // and U' respectively, 
which determine no in the Mott regions; ii) the regions 
of superfluidity between the Mott lobes shrink in size on 
the zJ/U-^/U phase diagram. If the coupling to the 
bath is strong enough, superfluidity is absent for small 
zJ/U regardless of it, illustrating that coupling to exter- 
nal degrees of freedom inhibits the formation of a phase 
coherent state. We discuss possible connections of our 
findings to measurements performed in the presence of a 
magnetic trap potential and suggest how the coupling to 
a bath might be detected. 

The starting point of our calculations is the standard 
BH Hamiltonian, and we assume that the bosons couple 
locally to the heat bath, so that the total Hamiltonian is 



H — Hbos + Hbath + H coup , 



(1) 



where H COU p describes coupling to the bath. The boson 
Hamiltonian Hb os consists of a local piece Hq and a hop- 
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ping term Hj: 

Hbos 



Hi 



bath 



-J^2(a\a,j + atoj) - fi fij 

+ J2^h i (n i -l)=Hj + H , (2) 

i 

^e a b\J)i a , (3) 

ioc 

Ys^^la+ka)^- (4) 



fii = ajaj is the on-site number operator, and the nota- 
tion indicates that the summation is restricted to 
nearest neighbors only. The operators b[ and bi a create 
and annihilate the bath degrees of freedom which have an 
energy spectrum e a . This form of coupling, in which the 
bath degrees of freedom are coupled to the on-site density 
■hi, is most suitable to take into account fluctuations of 
the chemical potential. g a characterizes the strength of 
coupling, which we assume to be purely local and Ohmic: 



rye Q exp{-£ Q /A}. 



(5) 



where ry is the strength of the coupling to the bath, and 
A is the energy cutoff for the bath degrees of freedom. 
The inverse cutoff 1 /A represents the time scale for iner- 
tial effects in the bath [23j. We note that the model we 
consider is similar to that of trapped impurity atoms in 
a Bosc-Einstcin condensate considered in Rcf. 



24|. 



We focus on the mean-field phase diagram that, in the 
absence of the bath, can be obtained by decoupling the 
tunneling term Hj and finding self-consistently the con- 
dition for a non-zero mean-field superfluid order parame- 
ter ^> b ■ Our main goal here is to elucidate the qualitative 
effects of the coupling to a bath, hence our use of mean 
field theory - future work will investigate effects beyond 
mean field theory. Identical results are obtained if one de- 
rives the partition function with the effective action that 
describes fluctuations of ^b close to the transition line 
in the Mott phase. Following the line of Refs. [kJ |25| |. 
we write down the effective mean-field Hamiltonian: 



Hmf = H 



H 



bath 



H 



coup 



(6) 



The optimum value for ^b = zJ(di) and is in fact pro- 
portional to the superfluid density, with z being the co- 
ordination number. This can be seen, if we add and 
subtract Hmf from H , so that the mean-field value of 
the ground state energy per site is [T3| 

E = E MF (^ B ) - zJ{a\){a t ) + (a t )V B + (4)*b. (7) 

Emf(^b) is determined by the occupation number 
no(/i/U) in the corresponding Mott region, but here it 
is also a function of the parameters r\ and A characteriz- 
ing the bath. Since terms involving ^fg and ^ B can be 



treated as a perturbation in the Mott-insulator phases, 
we can expand Emf in powers of ^b and ^ B l^Sl 



E MF = Ej® + X \*b\ 2 + 0(|*s| 4 ). 



(8) 



In general x is a function of all parameters characteriz- 
ing the local part of the Hamiltonian. Self-consistently 
minimizing Eqs. ([7]) and ([5]), we obtain that ^b has a 
positive non-zero value once the condition 1/zJ = \ 
is satisfied [27]. It is not difficult to see from the per- 
turbation series in Eq. ([7]) that — x coincides with the 
zero-frequency Fourier component of the on-site Green's 
function for bosons at T — 0, 



iQ(t-t') = (fai(t)al{t')}. 



(9) 



In Eq. ([H]), ai{t) are Heisenberg operators and the angular 
brackets mean that the average is performed over the 
ground state of the system involving the bosons and bath. 

To proceed, we single out the local part of the Hamilto- 
nian Hi = H Q + H bath + H coup in Eqs. ©-Q, and employ 
the well-known Lang-Firsov transformation to obtain the 
new local Hamiltonian 

H[ = SHiS~\ S = exp J J2 9 f™M a ~ h a ) \ ■ (10) 

K ia a ) 

After standard manipulations [28j |. we obtain 



U' 



where 



(11) 



V = 11-2^,-, M' = M + £-- (12) 



We assume weak coupling to the bath which implies that 
U' is positive. The transformation Eq. (fT0| means that 
the Green's function Eq. can be re-expressed in terms 
of the operators a\ — S&iS^ 1 , defined in the Heisen- 
berg representation with respect to H[. Subsequent cal- 
culations are tedious but similar to those considered in 
Ref. El. Hence, 



ig(t-t') = ig (t-t')F(t-t'), 



(13) 



with 



F(t-t') = exp J -Y,(g a /e a ) 2 (l 

L a 



ie a \t—t | 



(14) 



and 



*0b(t-O = (n o + l)e < l*-K*'- t >0(t-t / ) 

+n o e i €+(*- t ')0(i / -t), (15) 
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where £± are the energies of particle-hole excitations 

£+=!*'- U'(n - 1) > 0, £_ = n' - U'n Q < 0. (16) 

At T — 0, the ground state occupation for bosons 
in Mott-insulating phases is determined as no — 
noip'/U') = Integer[///Z7'] + 1, or n = 0, if p! < 0, 
where U' and /i' are defined in Eq. (fT2|) . Using Eq. (J5]), 
one can easily show that Q^ a — » J cfe) 



F(t-t') = (l + iA|i-t / |)~ , ) 



(17) 



so that the equation for the phase boundary separating 
the insulating and superfluid phases is recast as 



1 

z~J 



iQ(u = 0) 
dx 



o (1 + zA)" 



(n + 1)6 



(18) 



where a; is an auxiliary variable of integration. The inte- 
gral in Eq. (JT5J) can be further expressed as a linear com- 
bination of confluent hypergeometric functions, albeit we 
have found it more convenient to integrate numerically 
to obtain the T = phase diagram. 

The form of the coupling in Eq. ((5|) implies that 
U' = U - 2nA and u' = u + r)A. This leads to two 
new important parameters n and nA/U in the problem 
as a result of coupling to a bath. We only consider the 
situation in which nA/U <C 1, so that the suppression 
of the bare interaction U is not too strong. This means 
that 77 can be comparable to unity only if A/U <C 1. 
Such a case corresponds physically to a narrow band of 
low-energy bath degrees of freedom strongly interacting 
with the lattice bosons, and is different from the picture 
in which A ~ 0(U) but the coupling is weak, n <C 1. The 
first scenario is close in spirit to the model considered in 
Ref. 24| , although in both cases novel features appear on 
the mean- field phase diagram. 

Figures [1] and display mean field phase diagrams 
calculated from Eq. (fP8]) in the limits of A/U > 1 and 
A/!J < 1 respectively, for several values of 77, such that 
nA/U < 0.25 in all cases. Since n 1S determined by u' 
and U', a given occupation number corresponds to lower 
fj,/U in comparison to the case in which n = 0. 

An important feature to note is the difference in the 
functional behavior of the phase boundary near the 
points of integer fi'/U' for a) weak bath coupling, rj < 1; 
and b) strong bath coupling, 77 > 1. From Eq. (fT8|) 
it follows that for 77 > 1, the integral over x converges 
even if we put or equal to zero, in contrast with 
the case 77 < 1. We consider these two cases separately 
and find the corresponding behavior close to the degen- 
eracy points. If one moves along the upper branch of 
the lobe corresponding to the occupation number no to- 
wards the point (///(/') = no from below, |£_| — > and 
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FIG. 1: Phase diagrams in the zJ/U-fi/U plane for A/U = 
3.0, and n = 0.0, 0.02, and 0.05. The internal parts of all 
lobes correspond to Mott phases with occupation number 
n = Integer[///(7'] + 1. 
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FIG. 2: Phase diagrams for A/U = 0.1 and n = 0.5, 1.3, and 
2.0. The meaning of regions is as in the previous figure. For 
77 = 1.3 and 2.0, the horizontal lines between zj = and zJ- 
separating the Mott regions with no differing by one are not 
shown. 



£ + — > U. At the same time, approaching the same point 
from above following the lower branch of the no + 1 lobe, 
|£_| — > U whilst £ + — > 0. Hence: a) If 77 < 1, it is easy to 
find the corresponding asymptotics for zJ from Eq. (|18p 
for small |£_| or £+: close to the point (///[/') = 710, 
zJ oc for the lower branch, and zJ oc £+~' 7 for 

the upper one. This behavior suggests a slightly nar- 
rower region of superfluid phase between the adjacent 
Mott lobes than would occur for 77 = 0. b) For 77 > 1, 
there is never a superfluid at zJ — and a finite value 
of z J is required to achieve superfluidity for any value of 
the chemical potential /_t. Denoting as z J_ and zJ+ the 
limits obtained exactly at (///[/') = no for the lower and 



4 



upper branches respectively, we find from Eq. (118]) . that 



d.r 



- [(n + 1) + (n a + 1 ± l)e 



-xUl 



1 



zJ± Jo (1 + xA) 



These two limiting values are not equal to each other: 
z.J + > zJ_, and for A/U <C 1 differ from each other by 
0(A 2 /U 2 ). At the lowest order in A/U, zJ± w (n + 
l)(n — \)(A/U). The horizontal line between zJ + and 
z.J-, albeit short, separates the superfluid phase from the 
Mott insulator having occupation number equal to no. 
There is also a line between z J = and z J- (not shown 
in Fig. [2]) which separates Mott states with occupation no 
and no + 1. These two features of the mean- field diagram 
are new and arise as a result of strong interaction between 
the bosons and bath. 

It should be emphasised that whilst our results are 
mean field expressions for the BHM coupled to a bath at 
zero temperature, they form the basis for further investi- 
gations of effects beyond mean field theory and at finite 
temperatures. We expect that the addition of fluctua- 
tions beyond mean field theory will modify our expres- 
sions for the phase boundary between Mott insulating 
and superfluid phases, although the result that a finite 
value of J is required for superfluidity at all « will be 
robust. At low, but finite temperatures, there will be 
depletion of the condensate and, as a consequence, en- 
largement of non-superfluid parts of the phase diagram. 
Studies of the BHM (in the absence of a bath) indicate 
that the changes from the T = case are largest around 
the points where fJ>/U is an integer [13| . Hence, in anal- 
ogy with Ref. [l]| we expect that, for rj < 1, superfluidity 
will be suppressed around the points (fi'/W) = n . For 
n > 1, in its turn, the transition between no and no + 1 
Mott states at small J becomes a crossover, while the 
superfluid regions on the whole will be shifted towards 
larger zJ/U. 

In addition to the general features we have outlined 
above, it is also of interest to connect our results to 
experiments in optical lattices. These experiments are 
performed in the presence of confining magnetic traps, 
rather than the uniform system considered above. This 
means that the trap potential — V(r) must be added to 
the chemical potential u. This leads to the appearance 
of the so-called "wedding cake" structure 25|, [30( con- 
sisting of concentric regions of Mott phases separated by 
shells of superfluidity in which atoms are mobile. Our 
results indicate that coupling to the bath decreases the 
thickness of the shells, thus facilitating the crossover from 
superfluid to insulator upon increasing V / J. 

It should be emphasized that the model discussed here 
is independent of other mechanisms of dissipation, such 
as excitation of the phase slips in a moving condensate 
observed recently [2C 



2l|. If a heat bath plays a role 



strength. Unfortunately, however, for the model consid- 
ered here, it is difficult to determine from the currently 
available experimental data how strong interaction with 
the bath may be. However, we believe that this particle- 
conserving interaction is not unrealistic on timescales rel- 
evant to experiments, and might be detected by perform- 
ing lattice depth modulation experiments similar to those 
discussed in Ref. [3l|]. I n these experiments, the lattice 
potential is modulated at a fixed frequency for a short pe- 
riod of time. The modulation effectively translates into 
variation of J. During this process, energy is transferred 
to the atoms resulting in heating and the broadening of 
the central momentum peak, indicating an increase in 
the non-superfluid fraction of the atoms as observed af- 
ter ballistic expansion once the trap is switched off. The 
presence of coupling to a bath may mean that a system 
heated quickly by modulation will cool with time. Hence, 
varying the time between the end of the modulation and 
the free expansion may reveal a decrease of the out-of- 
condensate portion as a result of relaxation towards the 
initial equilibrium state giving a window to access the 
bath coupling. 

In conclusion, we have performed an analysis of the 
effects of coupling the BHM to an ohmic bath on the 
mean field phase diagram of the BHM. We find that the 
changes from the absence of a bath depend on both the 
strength of the coupling to the bath and the width of the 
bath spectrum. We believe our work should be relevant 
to any physical system that can be modelled with a BHM. 

We thank C. Chamon, J. McGuirk, and Kun Yang 
for discussions relating to choice of coupling to the bath 
and dissipation in cold atom systems. This work was 
supported by NSERC. 



in experiments on cold atoms in optical lattices, then it 
would be highly desirable to have a way to determine its 
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